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We report on a systematic finite size scaling analysis of in-plane penetration depth data taken 
on Bi2Sr2CaCu2 0g+{ epitaxially-grown films and single crystals, and of in-plane and out of plane 
data taken on YBa2Cu3C>6.7 aligned powder. It is shown that the tails in temperature dependence 
Oh' of the penetration depths, appearing around the transition temperature, are fully consistent with a 

finite size effect. This uncovers the granular nature of these cuprates, consisting of superconducting 
homogeneous domains of nanoscale extent, embedded in a non-superconducting matrix. 



Since the discovery of superconductivity in cuprates by Bednorz and Miiller [1] a tremendous amount of work has 
been devoted to their characterization. The issue of inhomogeneity and its characterization is essential for several 
Q ■ reasons, including: First, if inhomogeneity is an intrinsic property, a re-interpretation of experiments, measuring an 
O ' average of the electronic properties, is unavoidable. Second, inhomogeneity may point to a microscopic phase separa- 
. tion, i.e. superconducting grains, embedded in a non superconducting matrix. Third, there is neutron spectroscopic 
Qh| evidence for nanoscale cluster formation and percolative superconductivity in various cuprates [2,3]. Fourth, nanoscale 
spatial variations in the electronic characteristics have been observed in underdoped Bi2Sr2CaCu20s+<5 with scanning 
tunnelling microscopy (STM) [4-7]. They reveal a spatial segregation of the electronic structure into 3nm diameter 
d . superconducting domains in an electronically distinct background. On the contrary, a large degree of homogeneity 
has been observed by Renner and Fischer [8]. As STM is a surface probe, the relevance of these observations for 
bulk and thermodynamic properties has been clarified in terms of a finite size scaling of specific heat and London 
penetration depth data, revealing L a b ~ 50A and L c ps 68A for the length scale of the superconducting domains 
along the c-axis and in the ab-plane, respectively [9]. Fifth, superconducting domains with L a t w 300A were revealed 
by x-ray diffraction in oxygen doped La2Cu04 single crystal [10]. Sixth, in YBa2Cu3C>7_,5, MgB 2 , 2H-NbSe 2 and 
Nb77Zr23 considerably larger domains have been established. The magnetic field induced finite size effect revealed 
lower bounds ranging from L = 182 to 814A [11]. 

This paper addresses these issues by performing a detailed finite size scaling analysis of the tails in the measured 
temperature dependence of the penetration depths, appearing around the transition temperature. Specifically we 
analyze the data taken on Bi2Sr2CaCu20s+5 in the form of epitaxially-grown films [12] and high quality single crystals 
[13,14], and on magnetically aligned YBa2Cu30g.7 powder [15]. The paper is organized as follows. Next we sketch the 
finite size scaling theory adapted for the analysis of penetration depth data. Then we analyze the experimental data 
£f) ' and establish the consistency with a finite size effect, uncovering the granular nature of these cuprates, consisting of 
homogeneous superconducting domains with nanoscale extent, embedded in a non-superconducting matrix. 

Supposing that cuprate superconductors are granular, consisting of spatial superconducting domains, embedded in 
a non-superconducting matrix and with spatial extent L a , Lb and L c along the crystallographic a, b and c-axis. in 
this case the correlation length £j in direction i, increasing strongly when T c is approached cannot grow beyond Li. 
Consequently, for finite superconducting domains, the thermodynamic quantities, like the specific heat and penetration 
depth, are smooth functions of temperature. As a remnant of the singularity at T c these quantities exhibit a so called 
finite size effect [16], namely a maximum or an inflection point at T Pi , where £j (T Pi ) — Li. There is mounting 
experimental evidence that, for the accessible temperature ranges, the effective finite temperature critical behavior 
. of the cuprates is controlled by the critical point of uncharged superfluids (3D-XY) [12,17]. In this case there is the 
universal relationship 
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between the London penetration depth and the transverse correlation length £* in direction i [17]. As afore- 
mentioned, when the superconductor is inhomogeneous, consisting of superconducting grains with length scales Li, 
embedded in a non-superconducting matrix, the £f's do not diverge but are bounded by 

e$ < li i^j^k. (2) 

A characteristic feature of the resulting finite size effect is the occurrence of an inflection point at T Pk below T c , the 
transition temperature of the homogeneous system. Here 
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and Eq.(l) reduces to 
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In the homogeneous case 1/ (Ai (T) Xj (T)) decreases continuously with increasing temperature and vanishes at T c , 
while for superconducting domains, embedded in a non-superconducting matrix, it does not vanish and exhibits a 
turning point at T Pk < T c , so that 



Ai (T)Xj(T) 



/dT 



= 
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Since the experimental data for the temperature dependence of the penetration depths is available in the form X a b 
and A c only, we rewrite Eq.(4) as 
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to derive estimates for the diameter of the superconducting grains along the c-axis and parallel to the ab-plane. Noting 
that in the homogeneous system the transverse correlation lengths diverge as 
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the critical amplitudes and associated critical properties of the homogeneous system can also be derived from a finite 
size analysis. From Eqs.(6) and (7) we obtain the relations 



Cab (T P J - L c , J? ab (T Pab ) e* (T Pab ) = L ab . 
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Indeed, the transverse correlation lengths cannot grow beyond the corresponding limiting length scales. Thus given 
T Pc and L c , T Pab and L abl the critical amplitudes ^Q ah and £q c are readily derived. The anisotropy follows then from 
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As the system feels its finite size when the correlation length becomes of the order of the confining length, a physical 
quantity O adopts the scaling form [18] 
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The point is that the scaling function / depends only on the dimensionless ratio L/^ (t, L = oo) and it does not depend 
on microscopic details of the system. It does, however, depend on the observable O, the type of confining geometry 
and on the conditions imposed (or not, in the case of free boundaries) at the boundaries of the domains. Considering 
the London penetration depth, it reduces to 
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For t small and Lk — > oo, ±y — > oo and 
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while for t = and 7^ 
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so that in this limit 
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As expected, a sharp superconductor to normal state transition requires domains of infinite extent. Moreover at 
t p , y = 1 and d(Xoi\oj/ (Aj (T) Aj (T 1 ))) /dt = 0. Accordingly, there is an inflection point at T Pfe . Since the scaling 
function 17 (y) depends on the type of confining geometry and on the conditions imposed (or not, in the case of free 
boundaries) at the boundaries of the domains, this applies to the amplitude gok as well. 

We are now prepared to analyze the extended and systematic penetration depth data of Osborn et al. [12] derived 
from complex conductivity measurements on cpitaxially-grown E^S^CaC^Os+a films using a two-coil inductive 
technique at zero applied field. In Fig.l we displayed the temperature dependence of the real part of the complex 
superfluid density which is proportional to 1/A^ 6 for films A, B and C with different doping level. For comparison wc 
included the leading critical behavior in terms of 



with the parameters listed in Table I. Apparently the data is inconsistent with such a sharp transition. It clearly 
uncovers a rounded transition which occurs smoothly and with that a finite size effect at work. Indeed the extreme in 
d$t (p) jdT clearly reveals the existence of an inflection point at T p below the bulk T c , where the transverse correlation 
length £* fc attains the limiting length along the c-axis (Eq.(9)). Using the estimates for T p , derived from the location 
of the extremum in d^t(p) jdT values for L c are now readily calculated with the aid of Eq.(6) and the parameters 
listed in Table I. The results included in this table clearly point to the nanoscale nature of the inhomogeneities along 
the c-axis. Nevertheless, due to the small critical amplitude of the transverse correlation length ^Q oh , which follows 
from Eq. (3) and the corresponding parameters listed in Table I, and the fact that L c is considerably smaller than the 
film thickness d, the critical 3D-XY critical regime is attained. 
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FIG. 1. Re (p) oc 1/X 2 ab and dRe (p) /dT versus T for films A, B and C taken from Osborn et al. [12]. The solid lines indicate 
the leading critical behavior of a homogeneous bulk system according to Eq.(16) with the parameters listed in Table I. 
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0.010 


0.011 


0.024 
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Tabic I: Estimates entering and resulting from the finite size scaling analysis of the in-plane penetration depth data 
of Bi 2 Sr2CaCu208+<5 films (A,B,C) and single crystals (SCI and SC2). The films are respectively, A overdoped, B 
slightly overdoped and C underdoped. 
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To extend the analysis further we consider the finite size scaling function entering Eq.(12) in terms of 
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In Fig. 2 we displayed this finite size scaling function for films A, B and C. The collapse of the data on a single curve 
indicates that the inhomogeneities, also differing in their extent L c , have nearly the same shape and are subject to 
the same boundary conditions. Indeed close to t = the data approaches the expected behavior (Eq.(12)) 
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and goc — 0.6. Moreover rewriting Eq.(15) in the form 
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this relation provides a consistency test of the estimates for go c , ^ ab and L c . From Table I it is seen that there is 
satisfactory agreement. 




t / Abs( t ) 

FIG. 2. Finite size scaling function g versus tj \t Vc \ for the films A (0)> B (□) and C (A) derived from the data of Osborn 
etal. [12]. The insert displays the log plot and the solid lines are Eq.(18). 



Next we turn to the microwave surface impedance data for the in-plane penetration depth of a high-quality 
Bi 2 Sr2CaCu 2 08+5 single crystal. In Fig. 3 we displayed the data of Jacobs et al. [13]. In analogy to the film data 
shown in Fig.l there is clear evidence for a rounded transition, giving rise to an inflection point. Using Eq.(6) the 
finite size scaling analysis yields the estimates listed in Table I (SCI). In this context it should be kept in mind 
that there is still a considerable uncertainty in the absolute value of the zero temperature in-plane penetration depth 
A a 6 (0), the estimates ranging from 1800 to 2690 A [19]. For this reason we considered in Table I \ a b (0) = 1800 and 
2500 A, leading to L c ~ 68 and 132 A, respectively. In any case, due to the small critical amplitude of the transverse 
correlation length ^ ab , the critical 3D-XY regime is attainable in both cases and in agreement with Fig. 3. 
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FIG. 3. Microwave surface impedance data for \ 2 ab (T = 0) /\ 2 ab (T) and d (X 2 ab (T = 0)/ \ 2 ab (T))/dT versus T of a high-quality 
Bi 2 Sr 2 CaCu2 08+« single crystal taken from Jacobs et al. [13]. The solid line is \ 2 ab (T = 0) /\ 2 ab (T) = 1.25 (1 - T/T c ) 2/3 with 
T c — 87. 5K, the dash-dot line its derivative, indicating the behavior of the homogeneous bulk system, and the dashed line is 
the tangent at the inflection point, T Pc w 87K. 

To substantiate the occurrence of a finite size effect and to clarify whether or not the rather small L c value is 
associated with a different shape of the superconducting domains and different boundary conditions, we derived and 
display in Fig. 4 the finite size scaling function. The parameters entering this derivation are listed in Table I. Although 
the curve adopts the generic shape, in analogy to the film data shown in Fig. 3, there is an essential difference. Indeed, 
go c = 1.6, entering Eq.(18), differs substantially from the film value c/oc = 0.6. 




y = t/Abs(t p ) 

FIG. 4. Finite size scaling function g versus t/ \t Pc \ for the in-plane penetration depth data shown in Fig. 3, derived from the 
data of Jacobs et al. [13]. The insert shows the log plot and the solid lines are Eq.(18) with go c = 1.6. 

To provide some hints concerning the extrinsic or intrinsic nature of the inhomogencities, we consider next the 
data taken on a high quality optimally doped I^S^CaC^Os+s single crystal of different provenance. In Fig. 5 we 
displayed the data for the sample with T c = 91K, obtained with the single coil inductance method [14]. In analogy 
to the data shown in Figs.l and 3 there is clear evidence for a rounded transition, giving rise to an inflection point. 
Invoking the finite size scaling analysis outlined above and the parameters listed in Table I (SC2) we obtain with 
Eq.(6) for the limiting length along the c-axis the estimate L c ~ 80A. Although the data is not very dense around 
T c , the resulting scaling function shown in Fig. 6 is consistent with the expected behavior of the finite size scaling 
function. The solid lines are Eq.(18) with g 0c — 1.2. 
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FIG. 5. Single coil inductance data for \ 2 ab (T = 0) /\ 2 ab (T) and d (\ 2 ab (T = 0) /\ 2 ab (T)) /dT versus T of a high-quality 
Bi 2 Sr 2 CaCu2 08+« single crystal taken from Di Castro et al. [14]. The solid line is \ 2 ab (T = 0) /\ 2 ab (T) = 1.9 (1 - T/T c ) 2/3 with 
T c = 91K, the dashed line its derivative, indicating the behavior of the homogeneous bulk system. 
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FIG. 6. Finite size scaling function g versus t/ \t Pc | for the in-plane penetration depth data shown in Fig. 5. The insert shows 
the log plot and the solid lines are Eq.(18) with go c = 1-2. 



Summarizing the in-plane penetration depth data taken on three E^S^CaC^Os+a films and two single crystals, 
we uncovered in all these samples clear evidence for a finite size effect, stemming from the tail in the temperature 
dependence of the in-plane penetration depth. We have shown that the scaling behavior of this tail is fully consistent 
with a finite size effect, arising from the finite extent of the superconducting domains along the c-axis (see Figs. 2, 
4 and 6). Noting that in these five samples L c is of nanoscale and varies from 68 to 137 Aonly, it is suggestive to 
trace the occurrence of domains back to an intrinsic phenomenon. Furthermore, from Table I it is seen that the 
go c value of the films differs substantially from that of the single crystals. Since <?oc depends on the shape of the 
superconducting domains and the boundary conditions at their surface, it becomes clear that the inhomogeneities in 
the films differ in an essential manner from those in single crystals. Moreover, given the estimates for g 0c , £g ah and 
L c , QocS^ab/ L c is readily calculated. The resulting estimates provide an independent check of the consistency of the 
finite size scaling analysis. Indeed, according to Eq.(15), this quantity should equal (A 0a h (0) /\ a b {T c )) ■ In Table I 
we observe satisfactory consistency. 

To substantiate the generic occurrence of the finite size effect we turn to the data taken on magnetic field c-axis 
aligned YBa2Cu30g.7 powder, where both, the temperature dependence of the in-plane and out of plane pene- 
tration depths, have been measured [15]. In Fig. 7 we displayed respectively, (1/A a {, (T)) 2 , d (1 / \ a b (T)) 2 / dT ver- 
sus T and d(l/\ ab (T) /A c (T)) /dT versus T. The solid lines are respectively, (1/A afc (T)) 2 = (1/A 0af) ) 2 |i| 2/3 and 
l/\ a b (T) /A c (T) = 1/ (AoahAo c ) |i| 2 ^ 3 , while the dashed lines are the respective derivatives. These lines indicate the 
leading critical behavior of the homogeneous system. The corresponding values for T c and the critical amplitude 
1/ (AoabAoc) are listed in Table II, while the parameters of the finite size scaling analysis are summarized in Table III. 
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Here we used Eqs.(6) and (7) to obtain the estimates for L c and L a b, the diameters of the superconducting domains. 
While L c is comparable to the values found in Bi 2 Sr 2 CaCu208+(5 (see Table I), L a b turns out to be an order of 
magnitude larger. Its value L a b — 592Ais consistent with the lower bound L a b > 576A, derived from the magnetic 
field induced finite size effect on the specific heat of YBa2Cu306.6 [H]- 




T(K) T(K) 

FIG. 7. (a) {1/Xab (T)) 2 (O) and A {1/X a b {T)f /AT (□) versus T for YBa 2 Cu 3 06.7 derived from the data of Panagopoulos 
et al. [15]. The solid lines is (l/A a ;,(T)) 2 = (1/A 0a b) 2 |t| 2 ^ 3 and the dashed line its derivative, indicating the leading critical 
behavior of the homogeneous system. The corresponding values for T c and the critical amplitude (1/Aoab) 2 are listed in Table II; 
(b) l/A o6 (T) /A c (T) (O) and d(l/X ab (T) /A c (T)) /AT (□) versus T for YBa 2 Cu 3 6 .7 derived from the data of Panagopoulos 
et al. [15]. The solid line is l/A a6 (T) /A c (T) = 1/ (X 0ab Xo c ) |t| 2 ^ 3 and the dashed one its derivative, indicating the leading 
critical behavior of the homogeneous system. The corresponding values for T c and the critical amplitude 1/ (AootAoc) are listed 
in Table II, while the parameters resulting from the finite size scaling analysis are summarized in Table III. 
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Table II: Estimates for T c , (A 0a h) , (Ao a 6A 0c ) , A 0c , 7 = Ao c /A 0a &, A a 6 (T c ) and A a h (T c ) A c (T c ) for YBa2Cu306.7 
derived from the data shown in Figs. 7. 
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Table III: Finite size scaling estimates for YBa2Cu3C>6.7 derived from the data shown in Figs. 7. 
(( (a) \t\~ v {Xoab/Xab (T)) 2 versus t/ \t Pc \ and (b) \t\~ v ((A 0afc A 0c ) / (A ab (T) A c (T))) versus t/ |t Pa J derived from the 
data shown in Fig. 7 for YBa2Cu30@.7. )) 

To complete the finite size scaling analysis and the evidence for superconducting domains with diameters L c and L a b, 
we displayed in Fig. 8 the scaling functions (A ab/A ab (T)) 2 versus t/ \t Pc \ and \t\~ u ((X 0a b^0c) / (\b (T) X c (T))) 
versus t/ \t Pab \. Together with Eqs.(13) to (14) we observe the required agreement with the behavior of the respective 
finite size scaling function. The near coincidence g 0c ~ g a6 — 0.5 implies nearly identical shape and boundary condi- 
tion along the c-axis and the ab-plane. Moreover, g 0c ~ 0.5 is close to the value in the Bi 2 Sr 2 CaCu208+5 films (Table 

I). To provide a consistency check of the finite size scaling analysis we note that, given the estimates for g k, ij^oi^oj 

and L k , the quantity ffofc^/^oi^Oj/^fc m rca dily calculated and should be equal to (AojAcy) / (Aj (T c ) Xj (T c ))(Eq.(15)). 
Tables II and III reveal satisfactory agreement. 
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FIG. 8. Finite size scaling function, (a) \t\ v (\ 0ab /\ a b (T)) 2 versus t/ \t Pr \ and (b) \t\- v ({\ Qab \ c) / (A a6 (T) A c (T))) versus 
t/ \t Pab \ derived from the data shown in Fig.7 for YBa2Cu306.7 with the parameters listed in Table III. The insert shows the 
log plot. 

To summarize, we have shown that the tail in the temperature dependence of the in-plane and out of plane 
penetration depth around T c , as observed in the experimental data considered here, is fully consistent with a finite 
size effect, arising from homogeneous nanoscale superconducting domains with diameters L a b and L c . Clearly this 
finite size effect is not restricted to the penetration depth but should be visible in other thermodynamic properties, 
including the specific heat, as well. In the specific heat it leads to a rounding of the peak and its consistency with 
a finite size effect was established for the data taken on YBa 2 Cu 3 07_5 high quality single crystals [17,20]. In these 
samples the domain size was found to range from 300 to 400 A. Although the investigations of Gauzzi et al. [21] 
on YBa2Cu306.g films with reduced long-range structural order clearly reveal that the size of the domains depends 
strongly on the growth conditions, we established their nanoscale size and their thermodynamic relevance in a variety 
of samples. 

We thank K. D. Osborn et al. for providing their experimental data. 
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